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ABSTRACT 

A model for the dc plasma used in carbon nanotube growth is presented, and one- 
dimensional simulations of an acetylene/ammonia/argon system are performed. The 
effect of dc bias is illustrated by examining electron temperature, electron and ion 
densities, and neutral densities. Introducing a tungsten filament in the dc plasma, as in hot 
filament chemical vapor deposition with plasma assistance, shows negligible influence on 
the system characteristics. 
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I. INTRODUCTION 


Carbon nanotubes (CNTs) have been receiving much attention for their potential in field 
emitter devices, displays, sensors, electrodes and other related applications. Growth on 
patterned substrates for these applications is facilitated by chemical vapor deposition 
(CVD). Recently, plasma enhanced CVD (PECVD) has been investigated by several 
groups and found to yield well-aligned multiwalled carbon nanotubes (MWNTs) and 
multiwalled carbon nanofibers (MWNFs ). 1 ' 28 Whereas the MWNTs grown by PECVD 
appear like nanotube towers, the MWNFs are individual, freestanding structures, but are 
more defective with occasional bamboo-like or stacked cone inner core . 23 29 A variety of 
plasma sources such as microwave, " inductive, ' rf capacitive, and dc 
sources 13 " 16 ' 18 " 20 have been used to create a glow discharge of the feedgas mixture. A 
common substrate heating method is through a resistive heater. The use of hot filament in 
CVD, popular in diamond deposition, is also found in the CNT literature ; 1,2,4 but, a dc 
plasma through a pair of electrodes is normally added to take advantage of the electric 
field effect on alignment. Therefore, this approach is called dc plasma enhanced hot 
filament CVD (dc-HFCVD) or hot filament PECVD (HF-PECVD). 

With wide ranging applications for CNTs, techniques that allow controlled growth on 
patterned substrates as well as those for large scale production to meet possible future 
commodity market needs (for structural composites) will assume increased importance. 
While early reports focused on growth recipes and application potential, it is important 
now to develop an understanding of the growth processes in order to realize commercial 
applications. In this regard, modeling and simulation of the process is a valuable tool. For 
example, a volume-averaged model of an inductive plasma of methane and hydrogen 
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revealed 30 the extent of methane dissociation (95%) and the creation of higher order 
saturated hydrocarbons, such as C 2 H 2 , C 2 H 4 , C 2 H 6 etc., among several radicals and ions. 

3 1 32 

A kinetic model coupled with a two-dimensional reactive flow simulation analyzed 
the high pressure carbon monoxide (HiPCo) process for the bulk production of single- 
walled nanotubes and correlated process variables to the growth rate. 

In the present work, a dc plasma used in MWNT/MWNF growth is modeled and results 
from one-dimensional simulations are provided. The effects of dc bias as well as a 
tungsten filament on plasma and neutral properties are discussed. Since dc plasma 
sources are widely used now for CNT growth, the present exercise provides useful insight 
into the process mechanisms. 

II. MODEL 


The governing equations used to model the dc plasma reactor are as follows: 
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Equation (1) is a continuity equation written for electrons and each of the neutral and 
ionic species, p is the mass density of the species; subscript s stands for species number, 
and t is time. W s is the source term representing the production/consumption of species 
through reactions. W s j is a similar term for reactivity at the filament surface, though H 2 
dissociation is the only reaction considered on the filament. The last term on the right- 
hand side of Eq. (1) is a one-dimensional representation of species loss due to flow 
through the exit where subscript i denotes inlet and ris a residence time. J s is the species 
flux described by Eq. (5) and consists of a drift and a diffusion term. E is the electric field 
causing the drift (of electrons or ions); ju s is the mobility (of electrons and ions); D s is the 
diffusivity of species s in the mixture; M s is molecular weight of species 5; P is the 
system pressure; P s is the partial pressure of species s; and R g is the gas constant. 

Note that the plasma is not in a thermodynamic equilibrium and therefore, individual 
temperatures for electrons, ions and neutrals need to be used in the analysis. The ions are 
heavy and exchange energy efficiently with the gas, and their temperature is taken to be 
constant at 0.1 eV. Therefore, energy equations are solved only for the neutral gas 
[Eq. (2)] and the electrons [Eq. (3)]. T and T e are the gas and electron temperatures 
respectively. T e here simply represents 2/3 of electron mean energy and its use does not 
imply the assumption of a Maxwellian distribution. C v is the specific heat; the subscript n 
denotes neutral species; h is the enthalpy; at is the thermal conductivity; and AH r ° is the 
heat of reaction for reaction r. The terms on the left side of the energy equations account 
for convective and conductive transport of energy. The neutrals gain energy through 
elastic collisions with electrons ( Q c ), charge exchange collisions with ions (Qce) and 
from the filament {Of). The electrons gain energy from the electric field through Joule 
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heating ( q e J e E) and lose energy through elastic collisions with neutrals. q e in the Joule 
heating term represents the charge per mass of electrons, -e/m e . In both Eqs. (2) and (3), 
the last term represents the heat ot reaction from various inelastic collisions. 

The plasma is essentially quasi-neutral in the center of the discharge but there is a 
significant charge separation in the sheaths. Poisson’s equation [Eq. (4)] relates the 
charge and potential distribution throughout the plasma. T is the potential, So is the 
permittivity of free space, and p c is the charge density. 

Because of earlier interest in using HFCVD systems for the deposition of diamond films, 
beginning in the late eighties, different means of incorporating the effects of the filament 
into a gas model were developed. 33 ' 46 The filament has two main effects in the HFCVD 
system: (1) the conduction of heat to the gas and (2) the catalytic dissociation of H 2 on 
the filament surface. Much of the filament modeling 33 ' 41 involved assuming an 
experimentally measured filament temperature and either chemical equilibrium or near 
equilibrium conditions with a specified incubation time at the filament given its high 
temperature. The earliest works 33 " 37 did not account for the catalytic activity of the 
filament, but later research 38 41 then included filament catalysis either by assuming the 
degree of dissociation of H 2 to be given by the equilibrium or near-equilibrium 
dissociation fraction at the filament temperature and reactor pressure or by altering the 
incubation time to match experimental data. Further work 4 " 46 moved beyond the 
chemical equilibrium assumption and developed expressions for the rate of catalytic 
dissociation of H 2 on the filament but continued to specify a filament temperature from 
experimental measurements. In the present work, the filament is treated as a point 
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source/sink in the neutral energy and H 2 and H density equations. A separate filament 
energy balance equation, in contrast to the previous research, is fully-coupled to the gas 
equations and solved for the filament temperature: 

25380 "j A H° C (6) 

>/ )m Hi a 

The power produced by passing a current through the filament is dissipated by 
conduction to the gas, radiation to the walls and heat loss by H 2 dissociation on the 
filament surface, as given by the three terms on the right-hand side of Eq. ( 6 ), 
respectively, / is the filament current; R is the filament resistivity; A is the cross sectional 
area of the filament; C is circumference of the filament; a is the Stefan-Boltzmann 
constant (5.6705xl0‘ 8 W m ’ 2 K' 4 ); e is the emissivity of the filament; T w is the wall 
temperature; Z is the rate constant for H 2 dissociation on the filament surface; Pu 2 is 
partial pressure of H 2 ; and M/y, is the molecular w r eight of H 2 . 

Two separate conditions are considered for the filament: one in which the filament is 
clean and the second in which the filament is carburized. The resistivity of a clean 
filament as a function of temperature is taken to be : 47 

R = -2.3936 + 0.0257} + 2.476 x 10 ~ 6 7} 2 - 9.996 x 10”"7}\//fi • cm (7) 

and the tungsten emissivity 48 is given by e= 0.12 (7}/ 1000). Carburization is the process 
by which carbon species adsorb to the surface of the filament and are incorporated into its 
structure. The effect of this process is to increase the resistivity and emissivity of the 
filament and also poison the catalytic decomposition of H 2 at the surface. Zeiler et al 49 
show that within one hour their tungsten filaments were completely carburized to form 
W 2 C in a gas mixture of 1% CH 4 and 99% fb at a pressure of 22.5 Torr. Given that C 2 H 2 


= HT f - T)^ + sa(T 4 - T^+ZjP, 
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carburizes tungsten more readily than methane, 50 it is likely that filaments under our 
conditions will be carburized even more quickly. To model a carburized filament, the 
resistivity 51 of W 2 C is used: 

R = 80 [1 + 2 .586 x 10 ~ 4 (T f - 298 )], juQ. ■ cm ( 8 ) 

along with the W 2 C emissivity 52 given by £ = 0.185 (7/ / 1000). The effect of poisoning 
the filament is modeled by reducing the dissociation rate by a constant. 


The first term on the right-hand side of Eq. (6) represents the heat conduction to the gas 
and is equal to the Qf term in Eq. (2) which will only have a value at the filament 
location. Here, the heat transfer coefficient h depends on Reynolds number Re through. 


h = K g / Z) (0.43 + 0.48 A /Re o ), W / m 2 K ( 9 ) 

where K g is the thermal conductivity of the gas and D is the diameter of the filament. The 
last term in Eq. (6) represents the energy dissipated in atomization of H 2 on the surface of 
the tungsten filament. The rate expression for this process in pure H 2 was first 


investigated by Langmuir in the earlier 1900s and more recently validated. This 
expression would result in the following for IT// j in Eq. (1). 
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with W H /= - As with Q f , these terms are zero away from the filament location. 


The boundary conditions at the anode and cathode are as follows. Potential and gas 
temperature are specified at the electrodes. At the cathode, the electron temperature is set 
to 0.5 eV 55,56 and at the anode a diffusion cooling boundary condition 57 is applied, 
V/2 = eVHA Ions densities are extrapolated, and electron densities are determined by 
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assuming a zero density at the anode and a secondary electron emission coefficient at the 
cathode of 0.05. The neutral densities at the boundaries are determined by setting the 
mass flux to zero. 

The governing equations were solved using a code called SEMS (semiconductor 
equipment modeling software). This is a multi-dimensional, multispecies, 
multitemperature, chemically reacting flow code. 58 ’ 59 The governing equations are 
discretized using an implicit finite difference scheme. The resulting difference equations 
are solved using a Gauss-Seidel line relaxation method. The code was first validated by 
running a few dc argon simulations for which results are available in the literature. 60 
Then, both dc and dc-HFCVD cases here were simulated in one-dimension (axial 
direction perpendicular to the electrodes). Two-dimensional effects are more important to 
establish uniformity issues over large wafers. Given that now the research is academic 
and most growth experiments are done on 1-2 cm pieces of wafers, ID simulations are 
adequate initially to understand the plasma characteristics, which is the primary purpose 
here. 

An acetylene/ammonia feedgas system is studied here as it appears to be the most popular 
among the users of dc plasma and dc-HFCVD reactors. 18 neutral species (Ar, Ar + , H 2 , 
H, C 2 H 4 , C 2 H 3 , C 2 H 2 , C 2 H, C 2 , N 2 H 4 , N 2 H 3 , N 2 H 2 , N 2 H, N 2 , NH 3 , NH 2 , NH and N), five 
charged species ( Ar + , NH 3 + , NH 4 + , C 2 H 2 + , C 2 H 3 + ) and electrons are considered. A set of 
140 reactions is included in the modeling. 61 These consist of electron impact reactions 
(dissociation, ionization, dissociative ionization), charge exchange reactions, neutral 
reactions, and electron-ion recombination. Though the system can have conceivably more 
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species and reactions, the present set was derived after eliminating species of relatively 
negligible concentrations and very slow reactions. No surface reactions on catalyst 
particles leading to the actual growth are considered. These reactions and mechanisms are 
unknown. Nevertheless, the current exercise considering only the plasma phase is still 
valuable, given that the actual growth on a 1-2 cm piece of a wafer for 10 minutes only 
provides a very small sink for the active species. Therefore, the computed densities tor 
most purposes are within the accuracy limits of the various reaction rate constants and 
collision cross sections and should help to develop an understanding of the plasma 
process. 

The grounded anode and the biased cathode are separated by 1.5 cm. The substrate is 
normally placed on the cathode to receive the advantage of the large cathodic electric 
field. A 0.25 mm diameter tungsten filament is positioned 1 cm above the cathode for the 
dc-HFCVD simulations. The substrate temperature is taken to be independently 
controlled by a combination of a resistive heater and water cooling system. The substrate 
temperature is taken to be 700° C in all simulations. Other base conditions are as follows, 
pressure of 8 Torr, wall temperature of 450° C, and flow rates of 100 seem of argon, 
50 seem of NH 3 , and 25 seem of C 2 H 2 . The dc bias on the cathode is varied between 325 
and 750 V. The filament current is varied from 0 A (corresponding to pure dc discharge 
only) to 1 5 A using both clean and carburized conditions. 
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III. RESULTS AND DISCUSSION 


A. DC Discharge 

We discuss the results for pure dc discharges first. In all the figures, the cathode sits at 
the zero point axial location on the left of the plots, and the anode is located on the right 
at 1.5 cm. Figure 1 provides the relative abundance of various neutral species across the 
electrode gap with no filament and a dc bias at the cathode of -325V. As expected, Ar is 
the most abundant since its inlet mole fraction is 0.571 and it is an inert species. The 
concentration of the other two input gases C 2 H 2 and NH 3 are approximately 7.4><10 15 and 
4.3 X 10 13 cm ' 3 (volume-averaged) at the cathode compared to the input of 3.4><10 16 and 
1.7*10 ' 6 cm ’ 3 thus representing 78% and 99.7% dissociation, respectively. FL is 
produced in abundance through atomic hydrogen abstraction and recombination. Atomic 
hydrogen itself is produced from a variety of C 2 H 2 and NH 3 dissociation reactions. A 
stable hydrocarbon such as C2H4 is generated in the plasma and exhibits a density of 
~2*10 12 cm' 3 . It is interesting to note the complete dissociation of NH 3 even at this 
relatively low bias of -325 V (low for a dc discharge). Han et al 62 indicate that if NH 3 is 
replaced with N 2 , only carbon films result; pure C 2 H 2 without ammonia also did not 
provide any nanotubes. They speculated that the H supplied by NH 3 dissociation provides 
additional etching effects, thus enhancing the formation of nanotubes. Chhowalla et al} 9 
also emphasize the etching effect of NH 3 as an important factor in growing vertically 
oriented nano fibers. 

The effect of dc bias on electron temperature is shown in Fig. 2(a). The electron 
temperature at the edge of the cathode sheath is very high due to the energy gained from 
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the sheath electric field. The corresponding electric field profiles are also shown in 
Fig. 2(b). As the dc bias increases, the sheath field increases which leads to higher peaks 
for electron temperature. The high electron temperature will lead to high ionization rates 
and since ionization is endothermic, electrons lose energy following ionization. The drop 
in electron temperature, after the sheath and in the bulk, is due to this energy loss and 
heat conduction. The electron energy rises again in the anodic sheath due to local electric 
field. All this behavior is well known in the dc discharge literature. 57,60 The electron 
temperature in the bulk is about 1-1.5 eV depending on the bias. This compares with a 
Langmuir probe measurement of 1.5 eV reported in Ref. 19. The higher electric fields in 
the sheath [Fig. 2(b)] with increased dc bias are responsible for the better alignment of 
MWNFs as observed in Ref. 16 and 19. 

The effect of bias on ion and neutral densities is displayed in Fig. 3. The ion densities at 
the electrodes are the lowest since they recombine at the walls. The peaks seen outside 
the cathodic sheath (x = 0.1 - 0.15 cm) are due to the high electron temperature and 
density locally. In general, the shape of the ion density profile is dictated by the local 
balance among competing processes: production by ionization, loss by recombination in 
the bulk, wall loss and smoothing out by diffusion. Figure 3(a) illustrates that in the 
cathode sheath, regardless of the bias, the argon ion is always dominant followed by 
CiH 2 + and then NH 4 + . This phenomenon is consistent with the differences in the ion 
diffusivities, where Ar + has the lowest, NH 4 + has the highest, and C2H2 is in between, 
which results in differing rates of recombination at the cathode. The NH 4 ion is displayed 
here because it is the most dominant nitrogen-bearing ion given the rapid charge 
exchange reaction NH 3 + + NH3 — > NH 4 + + NH, (2.1xl(f 9 cm 3 /sec). 63 Figure 3(a) also 
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shows that at the lowest bias of 325 V the dominant ion in the plasma bulk is the NH 4 ion 
followed by C 2 H 2 + and Ar + . The dissociation of NH 3 at 750 V is higher than at 325 V and 
thus the NH 4 density decreases, and Ar + and C 2 H 2 + become the dominant ions in the 
bulk. The difference in dominance in the bulk and sheath at 325 V has to do with 
asymmetric charge exchange reactions given in Table I in which argon and hydrocarbon 
ions react with NH 3 to produce ammonium ions without a reverse mechanism. This effect 
is not noticeable at the highest bias because of the reduced ammonia density. 

In Fig. 3(b), the density of an incoming gas like C 2 H 2 decreases with an increase in dc 
bias as more of the feedstock dissociates. Atomic hydrogen product goes up following 
increased dissociation of feedstock as the bias is increased. Not surprisingly, the general 
trend is that the feedstock concentrations decrease and dissociated product concentrations 
increase when the applied dc bias increases. Interestingly, Chhowalla el a /. 19 observed 
that the actual deposition rate decreases with an increase in dc bias, possibly due to 
greater etching by the dissociation product, atomic hydrogen. 

Figure 4 displays the effect of bias on the gas temperature. At the lowest bias, the gas is 
mainly heated by gas-phase endothermic reactions that follow from the radicals that form 
from electron-impact reactions. At the highest bias, as show in Fig. 3(a), the plasma 
density is three orders of magnitude larger than at the lowest bias. Thus, ohmic heating 
from elastic collisions between neutrals and electrons becomes important as represented 
by the local maximum in the temperature profile near the cathode. 


12 



B. Filament Effects 


Next, we examine the effects of the filament on the system characteristics. Figure 5 
displays the effect of changing the filament current on the gas temperature at 325 V. The 
gas temperature is lower up to almost 2 A of current than for 0 A, but at 2 A and above, 
the gas temperature indeed goes up due to filament heating. It is important to recognize 
that the filament does not always heat the gas, given the multiple sources of heating 
(substrate heater, plasma heating) and the level of filament current itself Included in the 
figure is one carburization case at 15 A. The effect is as expected; for the same current, 
carburization produces a higher filament temperature resulting from the higher resistivity 
and poisoning of the filament. The higher resistivity increases the powei to the filament, 
and the poisoning reduces the energy lost through H2 atomization. The increased 
emissivity due to carburization does increase the amount of heat radiated but not enough 
to counter the increased resistivity or poisoning. Figure 6 shows the effect of the filament 
current on the densities of C2H2 and C2H4 along with atomic hydrogen; other species are 
not shown for the sake of clarity. Whether clean or carburized, the filament has no effect 
on any of the carbon-bearing species at the cathode where the sample for nanotube 
growth resides. The only effect it has at the location of interest is to increase the level of 
atomic hydrogen, at most by a factor of two. 

Figure 7 shows the electron density and gas temperature as a function of filament current 
at 750 V dc bias. At 15 A, the gas temperature in the bulk plasma is lower than that for 
zero current. This is because the gas actually heats the filament under these conditions. 
Here, the heating resulting from the intense plasma is larger than that of the filament. The 


13 



electron density simply is not affected by the filament current. This was also the case at 
the lower bias although not shown. 

Figure 8 shows the effect of the filament on the neutral densities at 725 V. The figure 
shows an interesting result when contrasted with Fig. 6. At the highest bias, the filament 
has little effect on any neutral species including atomic hydrogen. This can be explained 
by re-examining Figs. (5) and (7). At 325 V, the gas temperature resulting from filament 
heating is larger than the non-filament case everywhere in the reactor chamber. This is 
not the case for the 750 V bias where the gas temperature is only affected by the filament 
in the region of the filament. This is likely due to the role of ohmic heating near the 
cathode which is not affected by the filament and is substantially larger than the filament 
heating in the cathode region. Since there are no observed differences in atomic hydrogen 
density with filament current at 750 V, this suggests that the actual catalytic activity of 
the filament is not important in directly determining the cathode atomic hydrogen 
concentration except through its effect on the gas temperature. As stated previously, even 
in the case without significant ohmic heating at 325 V where the filament does affect the 
gas temperature near the substrate, there is little effect on the species of relevance to 
carbon nanotube growth. This perhaps explains the observation by Han et al . b 2 that the 
plasma intensity plays a more important role than the filament current in CNT growth and 
that well-aligned nanotubes can be grown without filament current. Hayashi et al , 64 also 
arrived at similar conclusions since they were able to grow well-aligned nanotubes on the 
entire substrate uniformly with a dc bias instead of just in the small area under the 
filament. 
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IV. CONCLUDING REMARKS 


Plasma enhanced CVD has emerged as a means to grow vertically aligned MWNTs and 
MWNFs. DC plasmas have been reported in the literature both with and without the aid 
of a hot filament. Here we have developed a model to analyze the plasma physics and 
chemistry as well as filament effects. One-dimensional simulations of an 
acetylene/ammonia system have been performed. The dc bias on the substrate plays an 
important role as it controls the sheath electric field affecting CNT orientation, extent of 
dissociation, and production of neutral species. Unlike CVD growth systems without 
plasma enhancement, the filament does not appear to have any significant effect if other 
mechanisms of regulating the substrate growth temperature are employed. We are 
currently working on adding surface chemistry, though rudimentary, and with this, 
additional parametric studies and 2D effects can be investigated. 
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Table I. Asymmetric charge exchange reactions 


Reaction 

Rate (cm 3 /sec) 

Reference 

C 2 H 2 + + H 2 C 2 H 3 + + H 

1 . 0 x 1 O ' 11 

Anicich 63 (1993) 

C 2 H 2 + + NH 3 — > NH 3 + C 2 H 2 

2. 14x1 O ' 9 

Anicich 63 (1993) 

c 2 h 2 + + nh 3 nh 4 + + c 2 h 

9.61xlO’ 10 

Anicich 63 (1993) 

C 2 H 3 + + H -> C 2 H 2 + + H 2 

6 . 8 xl 0 " n 

Anicich and McEwan 65 (199 7) 

C 2 H 3 + + NH 3 — > NH/ + C 2 H 2 

2.48xl0" 9 

Anicich 63 (1993) 

NH 3 + + H 2 -> NH 4 + + H 

4.4x1 O ' 13 

Anicich 63 (1993) 

nh 3 + + c 2 h 4 -» nh 4 + + c 2 h 3 

h— * 

4^ 

X 

q 

Anicich 63 (1993) 

nh 3 + + nh 3 -» nh 4 + + nh 2 

2.1 xlO ' 9 

Anicich 63 (1993) 

Ar + + C 2 H 2 -> C 2 H 2 + + Ar 

4.2xlO " 10 

Tsuji et al. 66 (1993) 

Ar + + C 2 H 4 -» C 2 H 2 + + H 2 + Ar 

2 . 2 xlO ’ 10 

Tsuji et al. 66 (1993) 

Ar + + C 2 H 4 -> C 2 H 3 + + H + Ar 

8.35xlO ’ 10 

Tsuji et al. 66 (1993) 

Ar + + NH 3 — » NH 3 + + Ar 

1.6x1 O ' 9 

Anicich 63 (1993) 
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Figure Captions 


1 . Concentration of various neutral species at -325 V cathode bias. Pressure of 
8 Torr; substrate temperature of 700° C; and 100/50/25 seem of AJ/C2H2/NH3. 
No filament heating is used. 

2. Effect of dc bias on (a) electron temperature and (b) electric field; conditions are 
as in Fig. 1. 

3. Effect of dc bias on (a) ion densities and (b) neutral densities for the same 
conditions as in Fig. 1. 

4. Effect of dc bias on gas temperature for the same conditions as in Fig. 1 . 

5. Effect of filament current on gas temperature at a dc bias of -325 V. 

6. Effect of filament current on neutral densities at a dc bias of -325 V. 

7. Effect of filament current on electron density and gas temperature at a dc bias of 

-750 V. 

8. Effect of filament current on neutral densities at a dc bias of-750 V. 
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Fig. 2(b) 
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